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Abstract 

A new approach is given for the implementation of boundary conditions used in solving the 
Mukhanov- Sasaki equation in the context of inflation. The familiar quantization procedure is 
reviewed, along with a discussion of where one might expect deviations from the standard approach 
to arise. The proposed method introduces a (model dependent) fitting function for the z" jz and 
a" la terms in the Mukhanov-Sasaki equation for scalar and tensor modes, as well as imposes the 
boundary conditions at a finite conformal time. As an example, we employ a fitting function, and 
compute the spectral index, along with its running, for a specific inflationary model which possesses 
background equations that are analytically solvable. The observational upper bound on the tensor 
to scalar ratio is used to constrain the parameters of the boundary conditions in the tensor sector 
as well. An overview on the generalization of this method is also discussed. 
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I. INTRODUCTION 



As is well known, the inflationary paradigm provides solutions to many cosmological 
problems such as the flatness problem, the horizon or causality problem, and also dilutes 
unwanted (and unobserved) relics jl-3]. It also provides a natural mechanism of producing 
primordial perturbations that seed the inhomogeneities of the universe % |5j- The basic idea 
is that the quantum fluctuations of a classically homogeneous scalar field, the inflaton, source 
quantum fluctuations of the spacetime metric (the inflaton will create density perturbations 
which will source the scalar fluctuations of the metric). During the process of inflation, 
this quantum fluctuation is amplified to become a classical fluctuation, and, at the end of 
inflation, the fluctuation in the metric induces the density fluctuations of matter that were 
produced during reheating . This primordial perturbation generated during inflation then 
is what gives rise to the formation of structure in the universe. 

In this story, the crucial quantities to be determined are the amplitude of the primor- 
dial density and tensor perturbations, as the growth of structure is dependent on their size. 
The subsequent evolution of the primordial perturbations can be inferred from careful ob- 
servations of the history of the growth of structure. From this we expect that the density 
perturbation produced by inflation to be (9(10~ 5 ). 

Therefore, it is important to be able to accurately compute these perturbations in order 
to either preserve or rule out a given inflationary model. For the usual models of inflation, 
Ricci scalar plus a single canonically normalized scalar field, there are two components that 
will determine this amplitude. As will be described in more detail below, the first is the 
form of a function, z"/z, which arises in the Mukhanov- Sasaki equation. This equation is 
satisfied by a mode function, Vk, which arises by a redefinition of the co-moving curvature 
perturbation in momentum space, TZk, upon having written the original action in terms of 
IZk. Since this equation is crucial to finding the curvature perturbation (the same equation 
is also obeyed by the tensor modes), it is essential that one accurately specifies its form. 

The second component is the input from vacuum selection, which is equivalent to a 
boundary condition for the Mukhanov- Sasaki equation. For example, one avenue of study 
has been to alter the initial state to lie away from the standard Bunch-Davies vacuum 6- 
Such alternative boundary conditions are typically chosen by conditions set at a given 
cut-off scale in either momentum or time. These choices will then manifest themselves in 
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physical observables (for example as new features in the power spectrum, or enhanced non- 
Gaussianity) which can allow one to gain knowledge of the initial state from observation. It 
should be mentioned that there exist arguments [15] that the Bunch-Davies vacuum might 
be the most probable vacuum to produce the correct power spectrum from the perspective 
of technical naturalness, although this does not eliminate the possibility of deviations from 
the Bunch-Davies vacuum. 

It is customary to characterize inflation as a period of time where the scale factor grew 
almost exponentially, a period called de Sitter or quasi-de Sitter inflation (exact exponential 
growth of the scale factor, a oc e Ht , with H constant, is technically de Sitter inflation, and 
nearly exponential growth is termed quasi-de Sitter). If the universe inflates as a power 
law manner a oc t p where r is the conformal time and p is a constant, than the solution 
of the Mukahanov- Sasaki equation is known. Notice that de Sitter inflation correspond to 
p — — 1. The solution for general p is given by a linear combination of Bessel functions. 
The calculation for the perturbation amplitude has been well established for such a case.jjij 
However, for most models of inflation, power law expansion happens only in a short period 
of time either at the beginning or the end of the inflation. Thus this commonly used 
approximation may not apply to all inflationary models. If one insists on using the equations 
derived from the power law limit, one runs the risk of possibly ruling out phenomeno logically 
viable models, or of preserving models that are ruled out by observational data. 

In the present work we would like to address the possibility that the function z" j z inside 
the Mukhanov equation deviates from the de Sitter limit, and how that may affect one's 
choice of boundary conditions. It may be that before some point t p that using the de Sitter 
limit is not consistent, and therefore placing boundary conditions at r v is more natural. 
We have thus expanded the standard method of computing the amplitude of the primordial 
perturbation for such a case. Our method applies to those inflationary models that do not 
behave with power law (at least partially) expansion with some specific constraints in the 
background evolution. 

The paper is organized as follows. In SecJTTl we review the standard calculation of pri- 
mordial density and tensor perturbations where the de Sitter limit is taken. We also explain 
the physical reasons behind the commonly chosen Bunch-Davies vacuum. In Sec JIIII we in- 
troduce a new method of vacuum selection by applying this method to a specific inflation 
model. The principles of generalizing this method to other models is also given. In Sec JIVI 
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we give our conclusions. 



II. THE STANDARD METHOD 



We will now outline the key ingredients for the calculation of the primordial perturbations 
by quantizing the comoving curvature perturbation as well as the tensor perturbation 
(for a recent textbook treatment and lecture notes see The theories we consider 

here contain an Einstein-Hilbert action and a canonically normalized scalar field, minimally 
coupled to gravity with an arbitrary self-interacting potential 



S = J d'x^rg jJ^fl - \g» u d,ad v a - V (er) j 
where k 2 = 8irG = 1/Mj> z and our metric signature is — h ++. 



(1) 



A. Scalar Perturbations 



We begin with the perturbed Friedmann-Robertson-Walker(FRW) metric including the 
most general perturbations 



ds 2 = a 2 (r){-(l + 2A)dr 2 - 2d i Bdx i dr + [(1 + 211) + didjH T ] dx'dx*} 



(2) 



Where A (r, x) , B (r, x) , 1Z (r, x) and (t, x) are small perturbations around homogeneous 
FRW metric. We will be concerned with calculating the scalar TZ, which is the gauge 
invariant comoving curvature perturbation. 

Variation of the action Eq.([T]) gives the Einstein equations and the scalar field equation 
of motion, which at the background level are 
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where a prime denotes a derivative with respect to the conformal time, r, while a dot will 
indicate a derivative with respect to the coordinate time, t. 
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Putting the solutions for the background evolution back into the Einstein-Hilbert action, 
and expanding the action to second order in the perturbations gives (setting k 2 = 1) 



S, 



(2) 



2 / d Xa Hi- 



ll 2 - a~ 2 (diliy 



(6) 



The above expression can be obtained using the gauge symmetry in the action to choose 
da = 0. One may define the Mukhanov variable 



v = zTZ, where z = a' 
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-2a 2 
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(7) 
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We have introduced the slow-roll parameter e. The second order action can be rewritten as 



(2) 



drd 3 x 



(v'y - {d lV f + ^ s 



(8) 



To quantize this action first define the canonical conjugate momentum of v, and then 
impose the usual commutation relation 

dL 



n, 



dv' 



v 



'■ [i ) (r ) x),n„(r,x')]=^ 3 »(x-x') 



(9) 



Henceforth we shall set h = 1. 

Now one performs a plane- wave expansion of the now quantum operator v(t, x) in Fourier 
space 

/(7 3 k 
^W{r)a^ + vl{r)aie^] (10) 

Requiring the canonical commutation relation between dk(r) and aj^r), [ak(r), ^/(t)] = 
(27r) 3 5^(k — k'), we will obtain the Wronskian condition for the mode function Vkij) 



l v k v k ~ v k v k) = -i 



The mode function in momentum space satisfies the Mukhanov- Sasaki equation 

v'l (r) + (k 2 - -) v k (r) = 



Upon introducing the second slow-roll parameter 



a 



Ha 



one can express — in terms of the first and second slow-roll parameters 



- = ^[l--r ] + e + -^--er ] + —e-—fj 



(11) 
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(13) 



(14) 



In general, Eq.( IT2|) with z" / z given in Eq.( IT4|) is difficult to solve analytically. For a 
special subset of general theories where e and 77 are approximately constants, the equation 
is analytically solvable. In this special case z" j z can be written as 

— = r-^ where v = — I — (15) 

z r 2 1-e 2 K ' 

and the analytic solution for (r) is given in terms of Bessel functions 

v k (r) = aV^Ju (kr) + (3y^Y u (kr) (16) 

Where a and /3 are two complex parameters. A well known example for this special case is 
that of power law inflation, a = ct p , where e = r/ = Pj y L . One then obtains z = \J 2c ^ +1 ^ -t p , 
which gives z" 2 or v = —p + 1/2, and the comoving horizon (aif) -1 = r/p. Pure de Sitter 
expansion is the specific case of p — — 1. For genuine de Sitter inflation, z vanishes, which 
eads to an exactly scale-invariant power spectrum that is now observationally disfavored 
1^ | . This implies that inflation must deviate from the pure de Sitter case. 

The solutions to Eq. (I12p can be written either as linear combinations of Bessel functions, 
J n (x) and Y n (x), or Hankel functions, Hn\x) and Hn\x) 

v k (r) = a^J_ p+ i (kr) + PV^Y_ p+ i (kr) (17) 
= (kr) + Py/fH^i (kr) (18) 



The Wronskian condition Eq. ffTT]) requires 



*a p* m 1 — 1 2 

a p — ap = or a — 

2 1 1 



2 



(19) 



When the solution is expressed in terms of Hankel functions, there is a natural place where 
the boundary condition may be imposed, that is when \kr\ ^> 1, or equivalently when the 
comoving wavelength is deep inside the comoving Hubble radius. The asymptotic forms of 
the Hankel functions become positive and negative frequency modes 

lim y/TH {2) AkT) = e~ ikT ; lim ^H {1) ,(kr) = e ikr . (20) 

The fact that in the far past the solution approaches those of Minkowski space can be seen 
in the behavior of z" j z displayed in Fig.([T]) 



6 



z'Vz 1.0- 




FIG. 1: A plot of p(p — 1)/t 2 for the cases: p = —2 given by the dashed line, p = — 1/2 given by 
the dotted line, and p = —1/3 given by the solid line. 

The vanishing behavior of z"/z in this asymptotic region ensures that solutions to Eq. (jl2p 
reduce to the Minkowski type in the far past. Thus, for these classes of inflationary mod- 
els there is a natural boundary condition that the solution should approach the positive 
frequency ougoing mode with no incoming modes 

„— ikr 



lim vJt) 



(21) 



This form is seen to match that of H^\[kr) in Eq.(l20l). The boundary condition Eq.( l2~TT) 
is known as the Bunch-Davies vacuum |6| . This has the effect of setting a = 1 and (3 = in 
Eq. ljTTl) . This appears as a natural choice, as one may think intuitively that at the beginning 
of time, all the particles (or positive frequency modes) should move forward in time, thus 

(2) 

eliminating the possibility of having a contribution from the H _ L (kr) term. 

P~ 2 

We would like to stress that although r — > — oo is a legitimate limit formally, but physi- 
cally there will exist a time where the physical wavelength will be comparable to the Planck 
length where quantum gravity effects should take place. This means in that region, the 
background evolution can no longer be treated classically. Due to the lack of a full quan- 



turn gravity theory, the boundary condition may be imposed at some later time where the 
physical wavelength is greater than the Planck length. The effect of setting the boundary 
condition at a finite time may be that the state does not reside in the ground state, but 
rather in some squeezed or distorted state ?], Isl, Q- 

For general energy contents of the universe, the form of the scale factor will no longer 
be a simple power-law (although during times where a single component is the dominant 
contributor to the stress-energy such as during matter or radiation domination, the power- 
law form is a good approximation). As an approximation, a standard analytic approach is 
to assume the expansion is approximately de Sitter, and therefore e « . Together 

with the smallness condition of the second slow-roll parameter r] = — <C 1 using Eq.( TT5"j) 
one finds 

i ( 22 ) 



z" a" 
z a 



Under this assumption, the solution of the Mukhanov- Sasaki equation is 

v k (r) = aVk^H^Jkr) + fiVk^H^Jkr) 



(23) 



where at and (3 are two complex parameters with four degrees of freedom, one of which is 
fixed by Eq.( fT9i) . two more by Eq. (1211) . leaving one irrelevant phase undetermined. With 
these conditions, the solution for the mode function is 

„— iki 



Vk (t) 



1 

KT 



(24) 



This leads to the well known relations for the scalar power-spectrum, P-ji, the spectral index 
n s , and the running of the spectral index a s 
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dhik 



k=aH 



Ak 3 e 



k=aH 



k=aH 

k d(H 2 /e) dr 



k 3 Pft dr dk 



(25) 
(26) 
(27) 



Note that each relation is to be evaluated at horizon crossing, k = aH (or equivalently 
t = r*), due to the fact that the perturbation is frozen when the comoving wavelength 
becomes stretched outside the comoving Hubble radius. 

The accuracy of this program strongly depends on the quality of how well the approxima- 
tion of the curve z" / z compares with the true function z"/z . Applying these equations to 
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models where z" / z is not well fit by the power law and slow-roll approximations may result 
in serious deviations from the observable predictions. This is precisely the problem we will 
address in Section IIII1 Before doing so we will first give an overview of quantization in the 
tensor sector. 



B. Tensor Perturbations 

The calculation for tensor perturbations mirrors that of the scalar perturbation. The 
starting point is once again the perturbed FRW metric with a transverse, traceless tensor 
perturbation, (we have omitted the scalar and vector perturbations seen previously which 
has no effect on the tensor mode evolution due to decoupling of the scalar, vector, and tensor 
sectors) 

ds 2 = a 2 (-dr 2 + (Sij + hij(x, r)) dx { dx j ) (28) 

Next, one expands the Einstein-Hilbert action to second order in the perturbation 

M 2 f 

S {2) = ^J drdx 3 a 2 [(^.(x, r)) 2 - (d^(x, r)) 2 ] (29) 



As in the scalar case, one expands hij in Fourier space in terms of plane waves with modes 

d 3 k 



K 



/d k \ 

where efj are the spin-two polarization tensors. One can then make the definition 



/4 (r) = \ah s k (r) (30) 



which leads to the action (here we have set Mpi = 1) 



iti) - k 2 - - (A* 



a 



k > 



(31) 



This action gives similar equations of motion as Eq.( [12l 

A*f (r) + (k 2 - -) Ail (r) = (32) 



a r 

For a scale factor of the power-law form, the calculation follows exactly as in the scalar case 
while the fit function in terms of slow-roll parameter is now 

— = - 9 4 where /i = —~ \- t:- (33) 

a t z 1 — e 2 



In the case of quasi de Sitter inflation, e ~ 0, the power spectrum for a single polarization 
of the tensor modes toward the end of inflation is 

2H 2 



|2 

I A*jfe 1 



a 2 



k 3 

k=aH 



(34) 

k=aH 



which differs from the form of the scalar result in that the slow-roll parameter e is absent 
in the denominator. The full power spectrum is then twice this (due to two polarization 
states) Ph = 2P t = 4H 2 /k 3 , which leads to a small tensor to scalar ratio r = P h j Pn = 16e 
when the slow-roll parameter is small. 

These results for the power spectra are obtained under the assumption that the expansion 
is de Sitter or very nearly de Sitter in the sense that Eq. (T22]) is true. To obtain a more 
accurate predication one must solve the Mukhanov equation on a model-by-model basis 
using the exact z" jz (or a" /a) numerically, along with choosing a proper boundary condition 
accordingly. In the next section we will institute such a procedure in order to quantize models 
where Eq. (|22j) is not a good approximation. 

III. VACUUM SELECTION 

The mode equation one needs to solve is 

< 2 - /(t)H = (35) 

One can obtain the analytic expression for /(r) from solving the background equations. For 
the scalar case /(r) is given by z"/z, while for the tensor case it is a" /a. In the de Sitter 
limit, its value is shown in Eq. (1221) . In general, the background evolution is not tractable 
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ytically, and for the few cases where the background solution is analytically solvable 



20h27I|. the r dependence in /(r) may be so complicated that finding an analytic solution 
for Eq. (l35l) becomes impossible. It is worth mentioning that the equation is equivalent to 
a time independent Schrodinger equation in one dimension when r is regarded as space 
coordinate and / (r) is the potential for the wave function (r) . This analogy will become 
apparent when we introduce our method of solving Eq. ( 135ft . which is nothing but the usual 
WKB approximation in quantum mechanics. 

Under the condition that the background evolution of the metric is known for a particular 



inflationary model, one can determine whether the approximation Eq. (j22j) will be applicable 
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for the model under consideration. If so, then the system may be solved using the standard 
approach outlined in the preceding section. However, this is not the case for a great number 
of inflationary models. One can understand that this is so because the mechanism used in 
stopping inflation may falsify the standard approximation. Additionally, the initial confor- 
mal time can not always be pushed back to negative infinity where one would impose the 
Bunch-Davies boundary condition as is the case when the expansion is truly power-law. 

The essential idea of our method is that, since the information one needs in order to 
calculate observables to compare with experiment is the value of the mode function at 
horizon crossing (which is much later than the asymptotic time r —> — oo), it may therefore 
be advantageous to place the physical boundary condition closer to the time where we require 
accurate information. Imposing the boundary condition at negative infinity may lead to a 
situation where the approximate form for /(r) in Eq. fl22|) has deviated greatly from the 
actual solution due to the lengthy intervening period of evolution. Although we are placing 
the boundary condition nearer the era of observable inflation, we will continue to mimic the 
idea of BD vacuum selection in that we look at the time r = p where the wavelength of the 
mode is deep inside the horizon, and the effect of cosmic expansion is relatively small. The 
situation can then be approximated as physics in Minkowski space. 

In this section we will first demonstrate the method with an explicit example before going 
on to comment on considerations on applying the method in general. 



A. A Specific Example 

The model we use as an example has background evolution that is analytically solvable 
[20I . This particular model provides an interesting picture where the Big Bang is connected 
to inflation with a specific time delay. The behavior of z"/z is very different from the slow 
roll plus de Sitter limit, while the term a" /a is asymptotically equal to the de Sitter limit. 
We will apply our method to obtain the scalar power spectrum and constrain the tensor 
power spectrum from observation. 

Our point here is to show that there exist examples such that z" / z can not be approxi- 
mated by the standard fit function, thus there is a need of introducing new method of solving 
Mukhanov- Sasaki equation. We would like to point out this model does not provide a mech- 
anism of stopping inflation, therefore, even though our method predicts the correct power 
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spectrum in a certain parameter space, the phenomenologically viable parameter space is 
expected to change when the stopping mechanism is introduced. Therefore the model in 
question may be viewed as a demonstration tool (due to the attractive property that it is 
analytically solvable at the background level) rather than a fully complete model. 
We begin with an action of the form Eq.([T]) with the scalar potential 

2 



Via) 



6 



K 



csinh 




b cosh 




(36) 



The potential contains the dimensionless free parameters b and c, and k 1 is the reduced 
Planck mass as before. The background solution for all possible combinations of c and b 



20] . We will use the case c = 646 > as an example to illustrate our 



has been classified in 
method. 

The behavior of the scale factor of this model is to initiate expansion at r = tbb ~ 0.92 
where the scale factor is exactly zero. Beginning at a later time, 77 ~ 2.87, there is an 
inflationary period, and finally, when r approaches r M « 7.4, the scale factor diverges. At 
this point the physical time, t = J a (r) dr, will also diverge. Of course the finite value 
Tbb ~ 0.92 is not physically significant since r can be translated by an arbitrary amount. 
The solutions for the background evolution of the scale factor a (r) and inflaton cr(r) are as 
follows 



a t 




cn 



1 + cn (±r) 



4 \ cn 



(37) 



and 
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(38) 



where r = 2 \cE\* r. E is a free parameter of this model which determines the scale of a, 
therefore E can be chosen so that a is normalized in the conventional way, a (today) = 1. 
In the present application we will choose E so that 



2 \cE\ 



1 or r 



T. 



(39) 



The term cn (|r) = cn {\t\^) is the Jacobi elliptic function[28|. Following the standard 
method, one needs to determine whether e and n are approximately constants. If so, following 
Eqs.( TL5]) and fl33|) . one can determine the fit function of z" jz and a" /a. Since the analytic 
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FIG. 2: The slow-roll parameter e as a function of conformal time. 

solution of this model is known, we can plot the exact curves for e and 77 as in FigF2] and 
FigS}. 

One can easily see that the asymptotic values of e and rj approach the constant values 
and 1 respectively. These values will result in the breakdown of Eq. (ll5|) for z"/z and 
2/(r — Too) 2 for a" /a. From the above expression we can also compute z" / z and a"/a. The 
result is shown in the following plots along with the fitting function 

The fit function plotted in Fig. (5) is the predicted fit function for a" /a which is 
2/(r — 7.4) 2 , while the fit function for Fig. (4) is a quadratic function 

f fit (r)=m(r-p) 2 + h (40) 

with parameters 

m = -1.2,p = 5.66, h= -0.23 (41) 

In the case of the scalar perturbation, since the formula for the standard fit function 
breaks down, we now introduce a new quadratic fit function according to the following 
rules. The parameter p is designed to match the point where /(r) has slope zero while h is 
designed to match the value of / (r) at r = p. The value of m is the only parameter in the fit 
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FIG. 3: The slow-roll parameter r\ as a function of conformal time. 
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FIG. 4: A plot of the function z"/z as a function of conformal time. The solid line is the actual 
value of the function z" /z, and the dashed line gives the fit function. 



14 



FIG. 5: A plot of the function a" /a as a function of conformal time. The solid line is the actual 
value of the function a" /a, and the dashed line gives the fit function. 

function, but it has to be adjusted such that the curve of //«(t) approaches the actual / (r) 
in the region of physical interest, namely the region where the mode exits the horizon. That 
is to say we want a better fit for r > p. The region r < p corresponds to the time before 
the beginning of observable inflation, or before about 60 e-folds before the end of inflation. 
The perturbation generated during the preceding periods have yet to re-enter our horizon, 
thus they have produced no observable effect. (Since the universe is currently accelerating 
again, it is not clear to what extent regions corresponding to time intervals before the last 
60 e-folds will contribute to observable effects in the future.) 

Once the fit function is chosen, one can then proceed to solve the Mukhanov-Sasaki 
equation with the new fit function. In the current example, the two independent solutions 
to Eq.f l35|) with / (r) = ff it (r) are the hypergeometric functions Ski and The general 
mode function is thus 

v k (r) = aSki (r) + 0S k2 (r) 
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where a and /3 are two complex parameters that can be parametrized by four real parameters 

a = r ie if \ /3 = r 2 e if?2 (42) 

written in the real parameters 

v k (r) = e iei (nS! + r 2 e i ^- ei )5 2 ) (43) 

r ii r 2-,d\i #2 are to be determined by the boundary conditions. It can be easily seen that when 
the quantity of interest is the expectation value of v k (t) the parameter 8\ is an irrelevant 
phase. 

The exact solutions Ski and S k 2 are 

S k i (r) = ei-^^hypergeom f^ -V^+Jf-h 1 ^ {t _ f\ (44) 

5 fc2 (r) = (r - p) e(-^ +p )^ T hyper geom (zl zh^LL^. ~ h *^(j- p f\ (45) 

V 4 V™ 2 / 

where the hypergeometric function is defined as 



7 / 7 \ -■ a a (a + 1) 2 . . 

hypergeom (a, 6, z) = 1 + -z + - - - - z + ... (46) 

The boundary condition Eq.f fTTl) will require a and (3 to satisfy 

z(a*/3-a/3*) if m < 



1 = i - v' k *v k ) = 

i(a*P -a(3*)e^ mp ' if m > 



or 

- = 2znr 2 sin (0 2 - 0x) = -iC (m, p) (47) 
Where C is a number which depends on the value of m and p. When m < 0, C — 1, 

/ — 2 

when m > 0, C = e~^ mp . In the specific example discussed in this paper m — — 1.2 < 
For completeness purposes, we have also included the expression for C for postive m. 

The next step is to choose a corresponding boundary condition. Since Fig. (4) shows 
z" jz < in the whole region of physical conformal time, 0.92 < r < 7.4, the solution to 
Eq.( l35l) will be a combination of oscillating waves. We introduce the WKB type solution 

v k (r)=A(r)e^ (48) 

where A (r) and <fi (r) are two real functions. Inserting this ansatz into Eq. (I35I) one obtains 

A" + 2iA'<p' + iA<p" - A (4>') 2 = -co 2 A (49) 
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where 00 (r) = \Jk 2 — / (r). This complex equation can be separated into two real equations 

A" - A ((f)') 2 = -co 2 A (50) 

(A^y = (si) 



The second equation can be solved easily 



^=4 (52) 



where q is an arbitrary constant. 

To solve Eq.f l50p we invoke the WKB approximation which assume A is slowly varying 
with conformal time, and then the second derivative term can be neglected. This approxi- 
mation is true when the potential, /(r), is slowly varying with conformal time. This is true 
when the function / (r) is at a local extrema, which corresponds to the very moment when 
we impose the boundary condition. Under this approximation the Eq. ([5T)l) can be solved 
analytically 

<f>(r) = ± fu (r) dr (53) 



The full solution can than be expressed as a linear combination of two modes 

gj J uj(t)cLt g— i J u)(t)cLt 

v k (r) = c 1 = hc 2 = — (54) 



In the vicinity of r = p, to (r) is a constant, cu* = y/k 2 — h. The boundary condition we 
impose is to choose the parameter a and /3 so that in the vicinity of r ps p, the solution is a 
linear combination of incoming and outgoing plane waves satisfying the usual Klein-Gordon 
normalization 

g— iuj*T gio;*r 

when r w p, Vk (r) ps a + b—= where, |a| — |6| =1 (55) 

The proposed boundary condition does invoke more parameters than the Bunch-Davies 
boundary condition since we are not setting a to 1 and b to 0, but instead keeping them as 
general parameters. As will be shown later, the phenomenologically viable parameter space 
requires a close to one and b close to zero. This means the boundary condition is close to a 
purely outgoing wave with smal amount of incoming wave, and thus close to the standard 
BD vacuum. It is interesting that a small amount of incoming waves is necessary to produce 
correct predictions. Unlike the Bunch-Davies boundary condition, which is imposed at the 
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beginning of time, our method imposes the boundary condition at a finite time, r = p, where 
there is no reason to assume there only exists an outgoing mode. One may criticize that our 
method introduces too many new parameters in the boundary condition. However, this may 
also be the case if z" jz can be approximated by Eq. (ll5p where the standard Bunch-Davies 
condition is meaningful. This is because transplanckian physics could significantly alter the 
initial condition such that the Bunch-Davies condition is not valid at all [7J • Thus, although 
more parameters are introduced in our method, the freedom is not necessarily more than 
the standard method. 



1. 



Now we express a and b in terms of real parameters a = a§e tda , b = b o e l0b , \a\ 
According to the proposed boundary condition, the mode function and its derivative are 
given by 



H 0)1 2 
K (p)\ 2 



-lujp 



2u 



-xuoae 



+ b 



lbjp 



-lujp 



iube iu3p 



1 

2 



bl + b oy T+¥ cos (a) 



— + OJ 



ft 2 , 



1 + bl cos (A 



(56) 
(57) 



where A = 2up — 9 a + 9b. Here we have assumed u = y/k 2 — h > 0. In the example we are 
discussing h = —0.23, thus this assumption is always true for this case. However, with a 
different fit function, it is conceivable that h could be positive, giving the condition k 2 > h 
in such an instance. The idea behind this choice is that at r = p, the equation becomes 
a harmonic oscillator with constant frequency u. For a small period of time near p, the 
solution should then approach the usual harmonic oscillator where the vacuum is given by 
the lowest energy state. From Eq. (l571) we find two more constraints that can be used to 
solve r 2 , and (62 — OA 



\ v k 



at 



p 

e 2 



K(p)\ 2 = Wf 



e 2 



1 

u 



1 + ^ cos ( A 



0J 



1 + bl cos ( A 



(58) 
(59) 
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The solutions can be categorized according to the sign of m. We find that if m < 



n 



\ 



1 b 2 + b 0V /l + b 2 cos (A\ 



2cu 



UJ 



' UJ 



j- + uj[b 2 ) -b \/l + b 2 ] cos(A 



0i) = - sin" 1 



1 + 4 (6§ + bt) sin 2 A 



(60) 
(61) 

(62) 



while for m > we obtain 



ri = e 2 p 



\ 



2w 



,,2 / / Cl? 



r 2 = e 



u \ bl — b \/ 1 + b 2 cos (A 



1 



sin 



l + 4(b 2 + b 4 )sm 2 A i 

Then for the mode functions we have for m < 
/ 



(63) 

(64) 
(65) 



v k (r) = 



b o+ b oV 1+b o cos(A) 



2w 



Si 



+e 



4| bg+b^J sin 2 A 



^+Lo[b 2 -b y/l + b 2 COs[A))S 2 



(66) 



/ 



and for m > 



/ 



tfc (r) = e^ 1 



2^ 



+e 



/ 

(67) 



With the above expressions, the mode function is uniquely determined (up to an irrelevant 
phase 9i) by the fit function ff it (r) along with the parameters b and (9 a — 9^) from the 
boundary condition. As usual, the mode function will be evaluated at horizon crossing, 
where k = aH. 

In order to produce an observationally consistent spectral index and its running, we find 
the parameters A and 6q are and 0.075 respectively . The mode function then satisfies the 
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FIG. 6: The scalar spectral index Tig SIS 3b function of conformal time. 




FIG. 7: The running of the scalar spectral index a s as a function of conformal time, 
boundary condition 

at t « p. v k (r) « 1.0028-= + 0.075== (68) 

\/2w a/2^ 



One can derive the spectral index as well as the running of the spectral index using definition 
Eq. (JSD as function of conformal time 

At t 7.05 the spectral index is n ~ 0.95, and the running of the spectral index is a ~ 0. 
This is within the allowed value of WMAP constraints 



Now, consider the tensor perturbation. The mismatch of the standard fit function is 
shown in Fig. (5). The fit function 2/(r — 7.4) 2 is a good approximation in the asymptotic 
region. This is because the first slow-roll parameter is approaching zero when r > 6. Despite 
the success in the asymptotic region, the validity of Bunch-Davies boundary condition is 
questionable. The reason is that Bunch-Davies vacuum imposes a condition at a fictitious 
conformal time r — > — oo, while in the physical universe actually starts its expansion at a 
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finite conformal time r ~ 0.92. The limit r — > — oo is inapplicable in our model thus the 
mode function evaluated from the evolution of such ,an ill defined boundary condition is not 
reliable. Thus, we shall keep the standard fit function but abandon the boundary condition 
that is imposed at an unphysical time. We will proceed by obtaining a general solution 
of the mode function, and then subsequently use the observational constraints to restrict 
the possible parameters. Since the power spectrum of tensor perturbations has yet to be 
observed, experimentally what we have to constrain is the ratio of tensor perturbation to 
the scalar perturbation. 

The most general mode function for the tensor modes is 



fi k (r) = aVkrH^hkr) + bVk^H^hkr) (69) 



l 3/2\ nl I ~ u v "3/2 

, -ikr / j \ e lkT 



1 - — + 1 + 7- (70) 

V2k V krj J2k\ kr 



From the definition fl3U) one obtains 



Ph=A ^ + 8 ^ { b ° ~ 6 °\A+^ cos ( A )) ( 71 ) 

Together with computed from Eq.f l68|) as well as the observational restriction on r, one 
can deduce an upper bound for ^26q — 2b a/1 + cos ^^jj • From the WMAP7 data 
r < 0.2, which corresponds to (2bl - 26 a/1 + b 2 cos (aJJ < 0.16. 



B. The General Case 



We would like to comment on how to generalize this program, such that it may be 
implemented for a given background evolution of an inflationary model. There is no reason, 
a priori, that either Eq. fll5p in the standard method, or a quadratic function, as introduced 
in the previous section, should be a good fit for z" / z in general. 

Our proposal in this paper has been to point out that one should not apply the standard 
fit function to all inflationary models without carefully examining whether it is actually 
a good fit. If the analytic solution for the Mukahanov- Sasaki equation is not attainable, 
one may wish to find a fit function that is analytically solvable. There exist a number of 
analytically solvable fit functions for z"/z, the quadratic function introduced in the previous 
section is just one simple example. 
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Another more sophisticated example is the quartic function which renders Eq.( |T2|) solv- 
able by linear combinations of the Heun triconfluent function. Unfortunately in practice, the 
introduction of a more complicated fitting function, such as the quartic function, may reduce 
the predictivity of the model as there will tend to be more parameters in the fit function 
that require matching rules to pin down their values. However the rule of thumb is simple, 
the employed fit function should resemble the actual curve at the region of interest, which is 
the moment of horizon crossing. After that, one should impose the boundary condition at a 
meaningful time. As we have demonstrated, using the WKB approximation at the moment 
when z" I z is flat is one sensible choice. Finally, one can use the current observational data 
for the tensor to scalar ratio to constrain tensor perturbation. However, since the power 
spectrum of the tensor mode has not been observed, the constraints from the tensor sector 
are not overly restrictive. 



IV. CONCLUSIONS 



In order to make predictions testable by observations, inflation needs not only a model, 
but suitable boundary conditions. Some models of inflation do not seem to fall within the 
realm where the standard boundary conditions may be naturally applied. With this in 
mind, we have discussed the introduction of an alternative method which generalizes the 
standard approach of computing the scalar and tensor power spectra. It is suggests that for 
those models whose background is analytically solvable, one should re-examine their power 
spectrum using our method and find how does their spectral index compare with the results 
from standard method. 

In general, this procedure will introduce additional parameters into the model, thus 
allowing more accurate phenomenology, with the usual drawback that introduction of more 
parameters decreases predictivity. This new method is implemented on a model-by-model 
basis, hence the generic effects of this approach have yet to be determined. For the specific 
example discussed in this paper, we explored an inflationary model which has analytically 
solvable background dynamics. We introduced a quadratic fit (of course, other models 
may require more complicated fitting functions)for the function z" / z which appears in the 
Mukhanov- Sasaki equation for the scalar modes, and imposed boundary conditions at finite 
conformal time, r. It was found that near r = 7.05 the spectral index and its running, both 
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fall into a phenomenologically acceptable range. This calculation gives an example for the 
implementation of our approach, although the model in question is not fully realized in the 
sense that it lacks a proper accounting for the cessation of inflation in order to produce the 
requisite amount of e-fold expansion. 

The capacity for altering the calculation, and thus the values, of observables predicted 
by inflation via this new approach is clear. It may therefore be possible that models which 
were hitherto discarded may need to be re-investigated in the framework of this method. 
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